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Abstract 

It is known (E.L. Green (1997), O. Post (2003)) that for an arbitrary m € N one can construct a 
periodic non-compact Riemannian manifold M with at least m gaps in the spectrum of the corre- 
sponding Laplace-Beltrami operator -Am- In this work we want not only to produce a new type of 
periodic manifolds with spectral gaps but also to control the edges of these gaps. The main result 
of the paper is as follows: for arbitrary pairwise disjoint intervals {aj,Pj) c [0, oo), j = 1, . . . ,m 
(m € N), for an arbitrarily small 6 > and for an arbitrarily large L > we construct a periodic 
non-compact Riemannian manifold M with at least m gaps in the spectrum of the operator -A^, 
moreover the edges of the first m gaps belong to ^-neighbourhoods of the edges of the intervals 
{aj,j3j), while the remaining gaps (if any) are located outside the interval [0, L]. 
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Introduction 

In this paper we deal with non-compact periodic manifolds. The ?i-dimensional Riemannian 
manifold M is called periodic if there is a discrete finitely generated abelian group F acting iso- 
metrically, properly discontinuously and co-compactly on M. Roughly speaking M is glued from 
countably many copies of some compact manifold M (period cell) and each y e F maps M to one 
of these copies. 

Let M be an ?i-dimensional periodic Riemannian manifold. We denote by -Am the Laplace- 
Beltrami operator on M. It is known (see e.g. 11231 ) that the spectrum cr{-AM) of the operator 
-Am has band-gap sti^ucture, that is 

oo 

cr{-AM) = [jjk{M), (0.1) 
k=i 

where J'k(M) = [ak,bk] c [0, oo) are compact intervals called bands, ak,bk / oo, = 0. In 

general the bands may overlap. The open interval {a,fi) is called a gap if (ff,yS) n cr(-AM) = 
and a,yS € o-(-Am)- 

The existence of gaps in the spectrum is not guaranteed: for instance the spectrum of the 
operator -Ar« = - d - in S." coincides with [0, oo). It is easy to see (cf. H) that in 1- 
dimensional case any periodic Laplace-Beltrami operator has no gaps. However in the case n>2 
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we have essentially another situation. Namely, E. B. Davies and E. M. Harrell II 161 considered 
the manifold M = W (n > 2) with a periodic conformally flat metric gij = aSij, where a = a{x) is 
a periodic strictly positive smooth function. The authors proved that a{x) can be chosen in such a 
way that at least one gap in the spectrum of the operator -Am exists. 

Further, E. L. Green |[T2l for any m € N constructed a periodic conformally flat metric in M? 
such that the corresponding Laplace-Beltrami operator has at least m gaps in the spectrum. 

Manifolds of another type were studied by O. Post in [24 1, where the author considered two 
different constructions: first, he constructed a periodic manifold M"^ (e > is a small parameter) 
starting from countably many copies of a fixed compact manifold connected by small cylinders 
(the parameter e characterizes a size of the cylinders), in the second construction he started from 
a periodic manifold which further is conformally deformed (the parameter e characterizes sizes of 
domains where the metric is deformed). For any m e N the existence of m gaps is proved for s 
small enough. These results were generalized by F. Lledo and O. Post ETl to the case of periodic 
manifolds with non-abelian group F. 

Also P. Exner and O. Post [7] proved the existence of gaps for some graph-like manifolds, i.e. 
the manifolds which shrink with respect to an appropriate parameter to a graph. 

We remark that a similar problem (i.e. the existence of gaps in the spectrum) was studied in 
im [TTJ [lH |30l for periodic divergence type elliptic operators in R", in [13] for periodic magnetic 
Schrodinger operator, and in 121 [T^ for periodic Maxwell operator. In these works the gaps in the 
spectrum are the consequence of a high contrast in the coefficients. We refer to the overview |[T5l 
where these and other related questions are discussed in detail. 

In the present work we want not only to construct a new type of periodic Riemannian mani- 
folds with gaps in the spectrum of the Laplace-Beltrami operator but also be able to control the 
edges of these gaps. Namely the goal of the work is to solve the following problem: for an arbi- 
trary finite set of pairwise disjoint finite intervals on the positive semi-axis to construct a periodic 
Riemannian manifolds M with at least m gaps in the spectrum of -Am (here m is the number of the 
preassigned intervals), moreover the first m gaps have to be "close" to the preassigned intervals, 
and the remaining gaps (if any) have to be "close" to infinity. 

Let us formulate the main result of the paper. 

Theorem 0.1 (Main Theorem). Let (aj,Pj) c [0, oo) f j - 1 , . . . , m, m e'N)be arbitrary pairwise 
disjoint finite intervals. Let 6 > be an arbitrarily small number, L > be an arbitrarily large 
number Let n e N \ {1). 

Then there exists an n-dimensional periodic Riemannian manifold M, which can be con- 
structed in the explicit form, such that 



Remark 0.1. In 1987 Y. Colin de Verdiere obtained the following remarkable result H: for ar- 
bitrary numbers = Ai < A2 < ■ ■ ■ < Am (m e M) and n e N \ {1} there exists a w-dimensional 
compact Riemannian manifold M such that the first m eigenvalues of the corresponding Laplace- 
Beltrami operator -Am are exactly Our main theorem can be regarded as an analogue of 
this result for the case of non-compact periodic Riemannian manifolds. 




m < m' < 00 



(0.2) 




(0.3) 
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Remark 0.2. Obviously it is sufficient to prove Theorem 0.1 only for such intervals {aj,/3j) that 
are nonvoid and their closures are pairwise disjoint and belong to (0, oo). For definiteness we 
renumber the intervals in the increasing order, i.e. 



< a\, aj < Pj < aj+i, j = \,m - \, a,n < P,„ < 



(0.4) 



Proving Theorem 0.1 we suppose that the intervals (aj,Pj) satisfy (0.4i. 



The idea how to construct the manifold M comes from one of the directions in the theory of 
homogenization of PDE's (for classical problems of the homogenization theory we refer e.g. to 
the monographs ll22l l27l l29l ). This direction deals with problems of the following type. Let M'^ 
be a Riemannian manifold depending on a small parameter e: it consists of one or several copies 
of some fixed manifold (we call it "basic manifold") with many attached small surfaces whose 
number tends to infinity as £ — > 0. On some PDE (heat equation, wave equation. Maxwell 
equations etc.) is considered. The problem is to describe the behaviour of its solutions as e — > 0. 
More exactly the problem is to find the equation on the basic manifold (so-called "homogenized 
equation") whose solutions approximate the solutions of the pre-limit equation as e ^ 0. 

Firstly the problem of this type was studied by L.Boutet de Monvel and E.Ya. Khruslov in ||2l 
where the behaviour of the diffusion equation was investigated. The asymptotic behaviour of the 
spectrum of the Laplace-Beltrami operator was studied in li5[ [T7l420l . in these works only compact 
manifolds were considered. 

Let us describe briefly the construction of the manifold M solving our main problem. We 
denote by Q.^ (e is a small parameter) a non-compact domain which is obtained by removing 
I" a countable set of pairwise disjoint balls D*^. (/ € Z", j = 1, . . . , m). It is supposed that 



from '. 



Df^. + ei and Df, . cz nf, = {x e R" : < Xa < e, Va). We denote by J"? the radius of the ball 
D^.. Let B^. {i € Z", j - 1, ... ,m) be an «-dimensional surface (we call it "bubble") obtained by 
removing a small segment from the n-dimensional sphere of the radius b^j. Identifying the points 
of dD'?. and we glue the bubbles B'r. (i e Z", j = I, ... ,m) to the domain Q.'^ and obtain the 
n-dimensional manifold M^: 



M'^ - Q'^ U 



uu 

iieZ" ;=1 



Bf 



The manifold (for m = 2) is presented on the Figure [T] 




Figure 1 : The manifold M" (m = 2). The period cell is tinted in more dark colour. 
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We equip M*^ with the Riemannian metric g'^ which coincides with the flat Euclidean metric in 
Q.^ and coincides with the spherical metric on the bubbles B^^: 



The manifold M"^ is periodic, the set 

m 



Mf 



Ff U 



7=1 



where Ff = Dq \ 



7=1 



+ ei 



is a period cell (for any / € Z"). 

We set = djE^^ if n > 2 and d'j = exp j if « = 2, - 
are some positive constants which will be chosen later. 



- bjS. Here dj, bj {j = l,...,m) 



We prove (see Theorem 2.1 1 that the spectrum cri-AM'^) of the operator - Am*: has at least m 
gaps when e is small enough (i.e. when e is less than some eq). We denote by {cr'^j,ju'^j) (j = 
1, . . . , m) the first m gaps, by J'^ we denote the union of the remaining gaps (if any): 



o-(-Am=) = [0, oo) \ 



(0.5) 



Then 



yj=\,...,m: limo-^^cT;, lim/i^=jU; 



lim inf J'^ = ca 



(0.6) 
(0.7) 



where the numbers cry, jij depend in a special way on dj, bj and satisfy the conditions 



0<Cri, O-j < jUj < O-j+i, j = l,m- I, 0-m<l^m< 



The set [0, oo) \ |J (cry,//y)j coincides with the spectrum of some operator ^ acting in the Hilbert 
space H = L2(W) © L2iM.",Pjdx), where pj (j = 1, . . . , m) are some positive constant weights, 

7=1, m 

by dx we denote the density of the Lebesgue measure. 

Remark 0.3. In the case when Q."^ is obtained by removing a system of balls from some compact 
domain Q. and m = I (i.e. the removed balls are equivalent, the attached bubbles are also equiv- 
alent) the behaviour of the spectrum of the Laplace-Beltrami operator with Dirichlet boundary 
conditions on dM"^ = dQ. was studied in [17J, also it was studied in [19] for another size of the 
removed balls, namely ei^ «c «; £ if n > 2 and exp^-^^ <^ d'j ^ s Q/a > 0) if « = 2. 
The same manifolds were also considered in (3^ where the behaviour of attractors for semi-linear 
parabolic equations was investigated. 

It was proved in ifTTl that the spectrum of the operator -A^^ (here ^ means the Dirichlet 
boundary conditions) converges in the Hausdorff sense (see the definition at the beginning of 
Sectionjs]) to the spectrum of some self-adjoint operator Jl^ acting in the space L2(n)®L2(r2,pJx), 
where p > is some constant weight. The spectrum cr(Jl) of the operator has the form 

a-iJi'^) = {o-} U [Af'' : A: = 1,2, 3...) U [/if"^ : k=\,2,3. 
where cr > is a point of the essential spectrum, the nondecreasing sequences /if ' , /if '^ belong 



to the discrete spectrum, moreover Um /l^ 



cr, lim A' 

k-^oo ' 



oo and A?'^ > fi, where = cr + crp. 



Thus, (o",ju)no"(J?l^) = 0, and, therefore, for an arbitrarily small 6 > the interval {cr+S,fi-6) does 
not intersect with the spectrum of the operator -A^^ when e = s(6) is small enough. A similar 
result is valid for the Neumann Laplacian -A^^^: the spectrum of the corresponding limit operator 



consists of the point cr and two nondecreasing sequences /if'^ such that lim A 



k = ^' 



lim A, 



oo. Moreover A^''^ = yu. It is important that cr, p are independent of the shape of the 



domain Q. and the type of the boundary conditions. These facts suggest that in the case O = R" the 
spectrum cr(-AM^) has a gap when e is small enough, and this gap is close to the interval {cr,fi). 



The proof of Theorem 2.1 consists of three steps. Firstly we prove that the set [0, oo) \ 



^ U coincides with the spectrum cr(^) of the operator Jl. Then we make the main step: 

m 

we show that for an arbitrary L ^ [j |jUyj the set cr(- A^'^) n [0, L] converges in the Hausdorff sense 

to the set cr{Ji) n [0, L] as e ^ 0. Finally, we prove that within an arbitrary finite interval [0, L] 
the spectrum ctC-A^e) has at most m gaps when e is small enough. Together with the Hausdorff 



convergence this fact will imply the properties (0.5 l-(0.7 1 (see Proposition 3.1 at the beginning of 
Section|3). 

We note that the metric g"^ is continuous but piecewise-smooth. However one can approximate 
it by a smooth metric g^f^ that differs from g*^ only in a small p-neighbourhoods of dB^.. Moreover 
when p - p{s) is sufficiently small then the spectra of the operator -A(M'=,g*^p) and the operator 
-A^E have the same limit as e ^ (here -A^M-'^g'^P) is the Laplace-Beltrami operator on M'^ 



equipped with the metric g'^f). For precise statement see Remark 4.2 at the end of the paper. 

In order to omit cumbersome calculations further we will work with the metric g^. 

Now, let (5 > be arbitrarily small number, L > be arbitrarily large number. It follows from 
Theorem 2.1 that there is such small e = s{6, L) that the structure of the spectrum o-{-Am'') is as 
follows: o-i-^M") has m gaps whose edges are located in ^-neighbourhoods of the edges of some 
fixed intervals {o-j,fij) (j = 1, . . . , m) while the remaining gaps (if any) belong to (L, oo). So we 



set M = M'^, E = e{6, L). In order to continue the proof of Theorem 0.1 we have to prove that for 



arbitrary preassigned intervals (ay, Pj) satisfying (0.4 1 it is possible to choose such dj, bj that 



CT: 



(0.8) 



We will prove this fact and present the exact formulae for the constants dj, bj (see Theorem 4.1 1. 

The paper is organized as follows. In Section [T] we recall some definitions and facts from the 
spectral theory for the Laplace-Beltrami operator. In Section [2] we construct the manifold and 
formulate Theorem |2 . 1 1 describing the behaviour of o"(-Am'!) as e ^ 0. Theorem l27T] is proved in 



Section|3] And, finally, in Section|4]we present the formulae for the parameter dj, bj. 



1. Theoretical background 

In this section we present the definitions and some well-known results related to the Laplace- 
Beltrami operator and periodic manifolds. For more details on the Laplace-Beltrami operator see 
e.g. ||28], for more details on periodic manifolds we refer to Il23l . 

Let M be an «-dimensional Riemannian manifold with the metric g. By gap we denote the 
components of g in local coordinates (xi, . . . , x„). 
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As usual we denote by L2(M) the Hilbert space of square integrable (with respect to Rieman- 
nian measure) functions on M. The scalar product and norm are defined by 



(m, v)z,2(M) ^ J" uvdV, \\u\\l2(M)^ V("' ")i.2W 

M 

where dV - ^Jdetgdxi . . .dxn is the density of the Riemannian measure on M. 

By C°°(M) (resp. Cq{M)) we denote the space of smooth (resp. smooth and compactly 
supported) functions on M. 

If the manifold M (possibly non-compact) has an empty boundary then we define the Laplace- 
Beltrami operator -A^ on M in the following way. By ?7m[m>v] we denote the closure of the 
sesquilinear form t/mL", v] defined by the formula: 

77m[m,v] = (Vm, Vv)l2(M) = J {Vu,Vv)dV (1.1) 

M 

with dom(77M) = C^(M). Here (Vm, Vv) is the scalar product of the vectors Vm and Vv with respect 
to the metric g: in local coordinates (Vm, Vv) = / , where g"^ are the components 

of the tensor inverse to ga^. The form rj is densely defined, closed and positive (by the way 
dom(77M) = H^{M) = {u e L2(M) : Vm € L2(M))). Then there exists the unique self-adjoint and 
positive operator -Am associated with the form timIu, v], i.e. 

(-Amm, v)l2(M) - fjM[ii,v] for all u € dom(AM), v € dom(?7M) 

For a smooth function u the Laplace-Beltrami operator is given in local coordinates by the formula 



-Amu = - ^^^^Vd^l^l (1-2) 



If M is a compact manifold with a piecewise smooth boundary dM we define the Laplace- 
Beltrami operator with Neumann (resp. Dirichlet) boundary conditions -Aj^ (resp. -A^) as the 
operator associated with the sesquilinear form fj^ (resp. fj^) which is the closure of the form rj^ 



(resp. 7/^) defined by formula (1.1 1 and by the definitional domain dom(?7j^) - C'^{M) (resp. 
dom(77^) = C~(M)). 

The spectra of the operators -Aj^ and -A^ are purely discrete. We denote by |/l^(M)|^^j^ 

(resp. |/l^(M)|^^^) the sequence of eigenvalues of -Aj^ (resp. -A^) written in the increasing 
order and repeated according to their multiplicity. 

Now we present the concept of periodic Riemannian manifolds. 

We say that the group F acts on the manifold M if there is a map Y x M ^ M (denoted 
(y, x) y ■ x) such that Vyi , 72 £ F, Vx e M one has (yi * 72) • ^ = y\ • {ji ■ x), where * is the 
group operation, and Wx e M one has id ■ x = x, where id is the identity element of F. 

The Riemannian manifold M is called periodic (or more precisely T-periodic) if a discrete 
finitely generated abelian group F acts on M, moreover 



• F acts isometrically on M, i.e. Vy e F: y- is the isometrical map. 



F acts properly dis continuously on M, i.e. for each x e M there exists a neighbourhood Ux 
such that the sets y • Ux (y e F) are pairwise disjoint, 
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r acts co-compactly on M, i.e. the quotient space M/T is compact. 



A compact subset M c M is called a period sell if |J y • M = M and M is a closure of an open 

yer 

connected domain D such that Vy € F, y id : D D y • D = 0. 

For convenience throughout our work we will use the same notation y for the element y € T 
and the corresponding map y- : M ^ M. 

By F we denote the dual group of F, i.e. the group of homomorphism from F into We 
remark that if F is isomorphic to Z" (as for the manifold M^, which will be considered in the next 
section) then F is isomorphic to the n-dimensional torus T" = {0 = {6i, . . . ,6r) € C" : Wa \0a\ = I}- 

Let 6 €r. We define the Laplace-Beltrami operator with 6-periodic boundary conditions - 
in the following way. By C^(M) we denote the space of functions u € C°°(M) satisfying 



u{yx) = 9{y)u{x) 

for each x e dM and for each y € F such that yx e dM. Then we define the operator -A^ as the 
operator associated with the form fj^ which is the closure of the form 77^ defined by formula (1.1 
(with M instead of M) and by the definitional domain dom(77^) = C^(M). 

The operator -A^ has purely discrete spectrum. We denote by |/l^(M)|^^^ the sequence of 
eigenvalues of -A^ written in the increasing order and repeated according to their multiplicity. 

For any 9 et the following inequality holds: 

if(M)<4(M)<if(M) (1.3) 

It turns out that analysis of the spectrum cr(-AM) of the operator -Am on the periodic manifold 
M can be reduced to analysis of the spectra cr(-A^) of the operators cr(-A^), 9 et. Namely one 
has the following fundamental result. 

Theorem. Let M be Y-periodic manifold with a period cell M. Then 



(T{-^M) = {^Jkm (1.4) 

kehl 

where J^kiM) = |/1^(M) : € f |, e N are compact intervals. 

2. Construction of the manifold 

In this section we construct the manifold M*^ and describe the behaviour of the spectrum 
cr(-AM'5) of the Laplace-Beltrami operator -Am-^ as e ^ 0. 

Let |d^^. : / € Z", j = I, mj be the system of pairwise disjoint balls in M." (n > 2) depending 
on small parameter e > 0. We suppose that: 

1) the balls D^j, j = I, m belong to the cube = |x e 1." : < ;Cq, < s, Va}; 

2) Vj - \, m: Ks < dist {d'^^., U ||J D^jj, where the constant a: > is independent of e; 

3) V/ e Z", V / - 1, m: D^. = D%. + si. 
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By x^.. we denote the centre of D^., by d'^. we denote the radius of If., (the third condition above 
lies that the radius of D^y depends only on the index ]). 

We denote by B'^. the truncated ?i-dimensional sphere (we call it "bubble") of the radius h'^:. 

^7 7 

fify - \(Qx,ei,...,Qn) : 01 e [0,2;:), ft e [0,7r), - 2,...,?i- 1, ft e [0^,7r]) 



Here 0^ - arcsin 



, where b^j (j = 1, . . . , m) are positive numbers satisfying b^j > d^.. 



Let us introduce in QF the spherical coordinates (^i, . . . , ft, r) with the origin at ;c?.. Here r is 
the distance to x^^.. Identifying the points (^i, . . . , 9n-.i,d'j) € dD'^j and (^i, . . . , ft_i, 0p € dB'!j we 
glue the bubbles B^^. to the perforated domain Q.^ and obtain an n-dimensional manifold M^: 



uu«? 

(6R" i=\ 



(2.1) 



The manifold M*^ is presented on the Figure [T] By ^ we denote the points of M"^. If the point x 
belongs to Q!^ sometimes we will write x instead of x having in mind a corresponding point in W. 

Clearly can be covered by a system of charts and suitable local coordinates {xi,. . . , x„) i-> 
X € Af^ can be introduced. In particular in a small neighbourhood of dB^j we introduce them in 
the following way (below by Ufj we denote this neighbourhood): 



Xt} — 



Xk = Ok, li - 1, . . . ,?i - 1, 
\r - d'j, X = (ft, . . . , ft_i, r) e n Ufj, 

[-b'j (ft - 0^) , x = {9u..., On-i, ft) e . n Ufj. 



(2.2) 



(that is dB'^j = {{xi, ...,x„): Xn = 0))). 

We equip with the Riemannian metric that coincides with the flat Euclidean metric on 
Qf' and coincides with the spherical metric on the bubbles B^.y This last means that in the spherical 
coordinates (^i, . . . , ft) the components of the metric g'^ have the form 



sip = Sap (b'jf Y\ sin^ ^k, a,l3=l,i 



/:=(!■+ 1 



(for a = « we set n sin ft:=l). Here 5^/3 is the Kronecker delta. 

k=a+l 

The metric g^ is continuous and piecewise smooth: in the coordinates (xi,...,x„), which 



are introduced above in the neighbourhood of by formulae (2.2 1, the components 
glp{x\, . . . , x„) of the metric g"^ have the form: 




Xn > 0, 
Xn < 0, 



(2.3) 



Where gl^^^ = 6,p(xn + d'^ \\ sin^ ft, gi^^^ = 6ap(b'^ sin^ + 0^ ]"[ sin^ft (2.4) 



k=a+\ V y / k=a+\ 

It is clear that as x,, - (i.e. on dB^) the coefficients g^ „ lose smoothness. 

Remark that can be approximated by a smooth metric g^^ that differs from g^ only in a 
small p-neighbourhood of dB^j, moreover when p = p{s) is sufficiently small then the spectra 
cr(-AM£) and o"(-A(M=,g':p)) have the same limit as e ^ (for more precise statement see Remark 



4.2 1. However in order to omit cumbersome calculations furffier we will work wiffi ffie metric g^ 



Remark 2.1. It is easy to see that the manifold M'^ can be immersed into the space R""*"' via the 
following map : ^ c R"+i (below x e R", z € R, (x,z) e R"+^): 

- if X = X € Q'^ then F%x) = (x, 0), 

- if X = (^1 , . . . , 0„) € B'lj then f^ix) = (xi , . . . , x„, z), where 

n n 

xi - + n sin 6i, Xk - + b"^- cos Gk-\ W sin Bi (k = 2, n), z - b%cos ©'^ - cos 0„) 

Note: one should not confuse (xi, . . . , x„) with the local coordinates introduced above in a neigh- 
bourhood of dBf.. 



Thus, maps B'^j onto the surface 5f . which is obtained by removing from the sphere S?. = 



X - x^^.r + (z - cos ^ (ftp^j the segment |(x, z) £ S'!j : z < 0) 



The map is a local homeomorphism, i.e. for any x € M'^ there is a neighbourhood U{x) c 
such that is a homeomorphism (and even diffeomorphism if x ^ IJ 5Bf ). If the surfaces 

B^^. (/ e Z", j = l,...,m) are pairwise disjoint (e.g. if b'j < d'j + ke/2) then is a global 

homeomorphism. Furthermore F'^ is an isometric map: if g'^ is a metric on which is generated 
by the Euclidean metric in R"^^ then g"^ coincides with the pull-back (F^)*g'^. 

Let the group = Z" act on M'^ by the following rule (below by y'^, k e Z" we denote the 
elements of F'^): 

- if X = X e Q'^ then maps x into the point yi^x - x + ks e QF, 

- if X = {9\,. . . ,Gn) € B^^. then maps x into the point y^x e B,+ytj- with the same angle 
coordinates {9\,. . . , 9n). 

Obviously is F'^ -periodic Riemannian manifold. For an arbitrary / e Z" the set 



Mf = Ff U 



U 

7=1 



, where Ff - < x € O'^ : x - e/ € Dq \ 



7=1 



(2.5) 



is a period cell. 

We assume that the radii of the holes and bubbles are the following: 



djE'-2, 



exp 



n>2 
n = 2 



b", - bjs 

J ■' 

where dj, bj {j - 1, . . . , m) are some positive constants (we choose them later in Section |4|). 
We will use the following notations: 



(2.6) 
(2.7) 



KS 



G^j-R^jUB^j, 



\x - x,.^.| 



= jxeQ'= 

a>„ is the volume of «-dimensional unit sphere 
9 



According to the notations introduced above in Section |lj we denote by {^f (Gf^) the se- 
quence of the eigenvalues of the operator -A^,, which is the Laplace-Beltrami operator in G?. with 

Dirichlet boundary conditions on S^.. It is clear that |/l^(G^.)} depends only on the index / 
One can prove (see Lemma |3 .21 below) that 



V; ^ 1, . . . ,m : 



where 



lim^i(Gf.) = cr,- 

£^0 ■' 



n = 2 
n>2 



2 d" ~U)n-\ 



(2.8) 



0, ii(Gfp does not 



Note that in spite of the fact that the diameter of G^^. converges to zero as e 
blow up as e — > 0. This is due to a weak connection between B'^.. and 7?^^.. 

We assume that the coefficients dj and hj are such that o", + (Tj if / j. For definiteness we 
suppose that crj < o";+i, j = I, . . . ,n - I. 

We introduce the Hilbert space 

//-L2(R") ®_L2(M",Pjdx) 

j=l,m 

where by dx we denote the density of the Lebesgue measure, the constant weights pj, j = I,. . .,m 
are defined by the formula 



Pj = (bjfoJn 



(2.9) 



Since lim ( /M) - 0, then p,- = lim e (here by | • | we denote the Riemannian volume). 
And, finally, let us consider the following equation (with unknown /I € K.): 



r{A) ^ 1 + V = 

(T: — A 



(2.10) 



It is easy to obtain (see the proof of Theorem 2. 1 1 that this equation has exactly m roots fxj (j 
1, . . . , m), moreover one can renumber them in such a way that 



O-j <IJ.j < O-j+i, j =\,m- 1, a-,n < Um < 



By the way if m = 1 then = cr\ + cr\pi (cf. Remark 0.3 1. 

Now we are able to formulate the theorem describing the behaviour of o"(-AM'!)■ 
Theorem 2.1. The spectrum cr(-AM£) of the operator —I^m''- has the following structure when e is 
small enough (i.e. when s < sq): 



o-(-AmO = [0,oo) \ 



(2.11) 



Here J''^ is a union of some open finite intervals (possibly J'^ - 0) and 



< (t\, (Tj<Pj< (T^+i, 7 = l,m - 1, < < //^ < MJ'' 
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Moreover 



V7= 1,. 



, , m 



f— >n J E—>0 ■' 



e-»0 

lim inf = oo 



(2.12) 
(2.13) 



m \ 

The set [0, oo) \ ^ IJ {(Tj,iJ.j)j coincides with the spectrum cr{^) of the self-adjoint operator ^ 
which acts in H and is defined by the formulae 



Ui 
U2 



-Ar,.M + ^ CTjPjiu - Uj) 

cri{u\ - u) 
0-2(^2 - ") 

0"m("m - m) 

dom(J?l) = dom(A]R«) ®_ L2(]R",py(ix) 

7=l,m 



e dom(j?l) 



(2.14) 



(2.15) 



We prove this theorem in the next section. In the last section we present the formulae for dj, 



bj which will ensure the fulfilment of the equalities (0.8 1 



3. Proof of Theorem 113] 

Before we prove the result in full detail we will sketch the main ideas of the proof. 
At first (Subsection |3.1| ) we prove the equality 



(r{:R) = {Q,oo)\ 



(3.1) 



In the main part of the proof (Subsections |3.2p3] ) we show that 

m 

for an arbitrary L > 0, L ^ |J the set o'{-Am'^) H [0, L] converges in the Hausdorff sense to 

7=1 

the set o-{J{) n [0, L] as e ^ 0. 

Let us recall the definition of Hausdorff convergence. 



Deflnition 3.1. The set : 



. converges in the Hausdorff sense to the set S c W as s ^ if the 



following conditions (A) and (B) hold: 



if A'^ € and lim = A then A € 



for any A e S there exists A"^ € 



3^ such that lim A"^ = A 



(A) 
(B) 



Property dAll is verified in Subsection 3.2 property dBl) is verified in Subsection 3.3 



In the last part of the proof (Subsection 3.4 1 we show that within an arbitrary finite interval 
[0, L] the spectrum cr(-AM=) has at most m gaps when e is small enough (i.e. when e < eo)- This 



fact and the Hausdorff convergence of cr(-AMO 1^ [0, L] to o-{J{) n [0, L] = [0, -L] \ ( |J {o-j,jij) 



imply the properties (2.11 1-(2.13 1. Indeed one can easily prove the following simple proposition. 
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U (o'y'/^pj, S = [0,L\\\\J (aj,/?;)], where L < oo and 



< Qfi, Oj < /3j < ffy+i, j = l,m - I, Um < L 
m'^ < m 



Suppose that the set converges to the set S in the Hausdorjf sense as s ^ 0. 
Then m'^ = m when s becomes small (i.e. when e is less than some eq) and 



Wj = l,...,m : lima^ = Qf;, limyS^-j6; 

£—►0 ^ e^O ■' 



3.1. Structure of cr{J{) 



Leti e C \ U cr; . LetF 



/i 



J[U -AU ^ F. Then u 



CTjU + fj 



e im{Jl - il), i.e. there is U 



Ml 



e dom(J?l) satisfying 



J 



and 



(Tj-A 

m n 

-A^nu - Ar(A)u - / + V 



7=1 



(3.2) 



where TiA) is defined by (2.10). Equality (3.2 1 implies that 



A € crm \ U [cTj] ^ AT {A) € (t{-Ak^) = [0, oo) 



(3.3) 



At first we study the function AT{A) on the real line. It is easy to see that AT{A) is a strictly 
increasing function on the intervals (-oo, cri), (cr^, oo), (cr,-, cr,+i), j = 1, . . . , m- 1, lim AT{A) = 

A—>±oo 



Figure 2: The graph of the function ATiX) (for m = 3). The bold intervals are the components of cr{J{) 
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±oo, lim XF{A) - +00, furthermore there are the points f^j, j = I,. . . ,m, such that 

nfij) = 0, J - l,...,m-l 
(Tj < Hj < O-j+i, j =l,...,m o-m<iu,„<oo 
m \ ( m- 1 ^ 

7=1 J V;=i 

Let us consider the equation Af{A) = a, where a e [0, 00). One the one hand it is equivalent 

(m \-l 

to the equation I (o"; - Pm+ 1 W - 0, where P„,+ 1 is a polynomial of the degree m + 1 , and, 

therefore, in C this equation has at most m + 1 roots. On the other hand it is easy to see that on 
[0, 00) the equation A'FiA) = a has m + 1 roots (if a = then these roots are 0,yUi, . . . ,//,„)■ Hence 
we obtain that the set € C : AT (A) > 0) belongs to [0, 00). 

The graph of the function AT (A), /I € R is presented on the Figure |2] 

Thus, we conclude that A e o-{J{) \ IJ |cry| iff i e [0, cri) U I (J [juj, o-j+\)\ U [pm, °°)- Since 

the spectrum cr(^) is a closed set , then the points crj (j - 1, . . . , m) also belong to cr(J?l). Equahty 
( |3.1| ) is proved. 

J.2. Property Q of Hausdorjf convergence 

We present the proof for the case « > 3 only. For the case n - 2 the proof is repeated word- 
by-word with small modifications in some estimates. 

Let A^ € c7"(-Am^) n [0, L] and lim A"^ = A. Obviously A € [0, L], thus, we have to prove that A 



belongs to cr{Jl). If A € [j o"; then this statement follows from (3.1 1. Therefore, we can focus 



on the case A ^ (j |cr,}. 

Let us consider the sequence sn c e, where en - jj^ ^ - 1,2,3... For convenience we 
preserve the same notation s having in mind the sequence sm- 
We introduce the following cubes in R": 

□ = |x e R" : < x„ < 1, Va) 
= |x e R" : Eia < Xa < E{ia + 1), Vck) , i - {ii,..., in) e Z" 

Since e"^ € N, then □ = |J u^, where 

= j/ e Z" : < < (e"^ - 1),Vq' 
Also we introduce the following set in M^: 



where is defined by formulae ( 2.5 1. 

In Section [2] we concluded that is F'^ -periodic manifold, the set M? is a corresponding 
periodic cell. On the other hand since e"^ € N, then is also F-periodic manifold on which the 
group F = Z" acts by the following rule (below by yu, k € Z" we denote the elements of F): 
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- if X = X e then maps x into the point y^x = x + k e Q.'^, 

- if X = {6i, . . . ,d„) e Bf. then yt maps x into the point y^x € Bj+^te-i / with the same angle 

I J 'J 

coordinates {9i,. . . ,6„). 

The set is a period cell. The boundary of M'^ is independent of e: dM^ = {x e Q.^ : x e da). 

Roughly speaking if e N then M'^ is not only "e-periodic" manifold but also "1 -periodic" 
manifold. To prove property ([A]l of the HausdorfF convergence it is more convenient to look at M'^ 
as F-periodic manifold (and to work with period cell M'^) since in this case we are able to utilize 
some ideas and methods developed in ll2l[3l [T7l[T9ll20]| . 

By Mq, (a = 1, . . . , 2n) we denote the components of dM'^: 

Ma = [xeQf : Xq, = and < x^ < 1, V/3 a} if o- ^ 1, . . . , n 

Mq, = ^x e Q.^ : Xa-n = 1 and < Xfj < i= a - if a = n + \, ... ,ln 

The faces Mq, and Mq,+„ (a = 1 , . . . , n) are parallel to each other and 



7,^,Mq = Mq+„, a - 1, . . . where - (0, 0, . . . , 1, . . . , 0) (3.4) 

T 

a-\h place 

Also we denote by Mq the corresponding faces of du. 

Since € cr(-AM'!), then there exists 6^ et such that A'^ € cr(-A^j). Since F is isomorphic 
to Z", then the dual group f is isomorphic to T" = [9 = {9],,..., 9^) € C : Va I^qI = 1}. For 
convenience hereafter by ff^ we will understand a corresponding element (9^, . . . e T". 

We extract a subsequence (still denoted by e) such that 

9F 9 = {9i,...,9n)eT" 

Let € dom(A^^) be the eigenfunction corresponding to A^, i.e. -A^^m^ = A'^u'^, + 0. We 
normalize by the condition ||m^||l2(M'') - 1, then IIVw^Hl^^jyjE) - A^. 

In order to describe the behaviour of as e ^ we need some special operators. From now 
on by C we denote a generic constant independent of e. 

We denote 



- X e : x e □ \ 



and introduce an extension operator 11^ : H^(M.^) H^{n) such that for each u e H^{M'^): 

Wuix) = u{x) for X € (3.5) 

l|n'^M||//ic) < ciimII (3.6) 

It is known (see e.g. Il22l Chapter 4]) that such an operator exists. 

By {u)b we denote the average value of the function u over the domain B c M*^ {\B\ + 0), 

i.e. {u)b = — J" udV^, where dV^ is the density of the Riemannian measure on M^. The same 

B 

notation remains for B c R.". 

If S c M'^ is a (« - 1) -dimensional submanifold then g"^ induces on S the Riemannian metric 
and measure. We denote by dS the density of this measure. Again by {u}-z we denote the average 

value of the function u over E, i.e. (u)b - — I udS'^ (here |E| = f dS'^). 
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We introduce the operators ffi : L2(M'^) — > L2(n) (j = 1, . . . , m) by the formula: 



Recall that □ = IJ Using the Cauchy inequality and (2.7 1 we obtain 



(3.7) 



In view of (3.6 1, (3.7 1 the norms ||n'^M'^||//i(Q), lin^w'^llLjCn) U = l,...,m) are bounded uni- 
formly in e. Usmg the embedding theorem (see e.g. Il28l Chapter 4]) we obtain that the sub- 
sequence (still denoted by e), the functions u e //'(□), uj e L2(n), j = I,. . .,m exist such that 

H^u^ M weakly in //'(□) and strongly in L2(n), — > m; weakly in L2(n) 
Moreover due to the trace theorem (see e.g. |[28l Chapter 4]) IT'^m'^, u € L2{dn) and 



n^u"^ u strongly in L2{dn) 



(3.8) 



Since e dom(A^^), then in view of (3.4 1 



Therefore, 



u^{x + Co) = 9%u^{x), X € M„, a = I,. . . ,n 



u{x + ea) = 9au{x), X € Mq,, a = I,. . . ,n 



Thus, u € dom{r]^). Recall (see Section [T]) that 77^ is the sesquilinear form which generates the 



operator -A^^. 

We also need some auxiliary lemmas. 

Lemma 3.1. For any j - I,. . . ,m: 



= u 



L2(n) 



(3.9) 



lim£«y 

Proof. We denote R^. = \x e Q.^ : d'^. + ^ <\x- x^.\ < d'^. + One has the inequahties: 

< lin'^Mll^^cp - ^" l<n'"'>nf P < Cs^im'u'WMn^r i e i"' (3-iO) 



(3.11) 
(3.12) 



which are valid for any € H (□□), j = \,. . . ,m. Inequality ( 3. 10 1 is the Poincare inequality, the 



inequality ( |3.1 1| ) follows directly from lfT9l Lemma 2.1], and the inequality ( 3.12| ) can be proved 
in the same way as inequality (2.2) from ||20lQ 

Equality ( |3.9| ) follows directly from ( 3.10| )-( [3TT2 i. The lemma is proved. 



'in \2Q\ inequality 1 3.12 1 with dR-lj \ Sfj instead of S'^j was proved. For Sf- the proof is similar We remark that in 
the case n = 2 inequality 1 3. 12 > is valid with | In e\ instead of s^^". 
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Lemma 3.2. 0For j = I, . . . ,m: 



limif(Gf.) = (r,- 

£—►0 -' 



where crj is defined by formula {2.8). 



Proof. Let v^^. € dom(A^£) be the eigenfunction corresponding to /li(G^y) such that {v'?j)b';. = L 

Instead of calculating vf . in the exact form we construct a convenient approximation v^. for it. 
We introduce the notations: 



B^j = {x = ieu...,en)e B^j : e„ € [0^, k/2]} 



5f. = {x - (01, . . . , 0„) € Bf. : e„ = nil] = dB^. \ dB^j 
Let the function v'?j be the solution of the following boundary value problem: 



-AG.v^ = OinGf^ (3.13) 



%^-0, =1 (3.14) 



'J 



Here by -A^s we denote the operator which is defined by the operation ( 1.2 1 and the definitional 



domain dom(A(j£ ) = |m : m = v\qs , v € dom(AM'!)|- For convenience from now on we use the 
notation -A instead of -A^e . It is easy to see that the function v^. is smooth in . and B^., the 
limiting values of . in the domains R^. and B^. coincide on dB'^., the normal derivatives satisfy 
the condition "g^ + ^ - 0- 

Due to the symmetry of one can easily calculate v'Tj (recall that we consider the case n >3): 

vUx) = \ ' " ' . (3.15) 

e„ 

where F(0„) = ^(sin'~" if/)di(f and the constants A^, B^, are defined by the formulae 

n/2 

(j'r)""^ A'^ A'^ 
A^ = B^- q = («-2)— ^ (3.16) 

l-(^) -(«-2)F(0p(^) (t) 
We redefine v?. by 1 in S^^. \ . preserving the same notation. 



Direct calculations lead to the following asymptotics as e ^ 0: 



where cTj, pj are defined by formulae (2.8 1, (2.9 1. 



^This result was given in 1171 without a justification. In the current work we present a complete proof. 
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We define the function v,- e dom(A^s ) by tlie formula 



vf 



I J 



1 + (vf .(x) - 1)0) , X = (01 , . . . , 0„) € Bf . 



(3.18) 



1, 



Here <1>(0„) is a twice continuously differentiable non-negative function on [0, oo) equal to 1 as 
< < njA and equal to as 9n > n/2. We have the following asymptotics as e — > 0: 



lim^-"(llvf,-l|li,(,.) + ||vf,|li,(«..Uo 
It follows from the min-max principle (see e.g. 1261 ) that 



(3.19) 
(3.20) 



^i(Gf.) = . ^ < 



(3.21) 



Note, that this automatically gives the inequality lim A\ (Gf ) < cry. 



We present the eigenfunction vf . in the form 



e-»0 



V . = V- ■ + W- ■ 

ij 



(3.22) 



Let us estimate the remainder wf^.. One has the following estimates for the eigenfunction 
the proof see {i^ Lemma 4.2]): 



vf . (for 



l|vf,|P 



< Ce' 



,n+2 



(3.23) 
(3.24) 



Using ( 3.20 1, (3.23 1, ( 3.24 1 we obtain 



^""ll<illL(G^.) ^ 2£-" (l|vf..||?,,,., + ||V^ "2 



/yM,,(«.) . M^,yM,,(«.) + l|l-vf^||i^(5.) + ||vf^-lii2 



(3.25) 



Substituting (3.22 1 into (3.21 1 and integrating by parts we get 



(Avfy,wOL2(Gfp 



+ ||VV 



e Ii2 



l|vf,-llL2(Gfp 



- 1 



(3.26) 



Taking into account ( 3. 17 1, ( 3. 19 1, ( 3.23 1, ( 3.25 1 we conclude that ( 3.26 1 implies 



e~"||Vw 



e ||2 







(3.27) 



It follows from pTTl), (|37T9ll, KISI), (|3:27ll that lim (Gf ,.) = o-j. The lemma is proved. 

1 11 11 r 1 t e^O ' 



Lemma 3.3. For j = I,. . . ,m: 



lim (Gf.) = oo 

£-»0 -' 



17 



Proof. Let G*! be an ?i-dimensional surface embedded into R"'*'^ (below x eR", z e R): 



G^ - U 



where 

R^ - {(x,z) e ]R"+^ : £~^J^ < Ui < k/2, z = 
- Ux,z) € R"+i : + (z - Z^^- cos 0^)^ = {bjf, z. > 



We equip G^ with the Riemannian metric induced by the Euclidean metric in R"+^ By dV we 
denote the density of the Riemannian measure on G*;. Thus, is the £~'-homothetic image of 

Evidently one has the following relation between the spectra of -A^^ and -A^^: 

u J 

\fken: if (Gf,.) = e-'if (Gfp (3.28) 

We denote 

R - {(x,z) € R"+i : \x\ < k/2, z = o) , By - ((x,z) e R"+i : \xf + z^ = (bjf] 
Further we will prove that 

V;t e N : if(G^) ^ Ak (3.29) 

where {/l/t)i:gN the eigenvalues of the operator which acts in the space L2(R) ffi ^2(B;) and 
is defined by the formula 

^ \0 Ab, 

Here the eigenvalues are renumbered in the increasing order and with account of their multiplicity. 
One has Ai = Ai(Bj) = 0, A2 = min{/lf(R),/l2(By)} > 0. Therefore, in view of (3.28 1-(3.29 1 

lim /l?(G^.) = 00. Thus, to complete the proof of the lemma we have to prove ( 3.29 1. For that we 

£^0 'J I 

use the abstract scheme proposed in the work ITSl . 

Theorem [16\. Let "K^, be separable Hilbert spaces, let Ji"" -.n-C ^ "K^, : -K" ^ be 
linear continuous operators, imj^^ CV c "H^, where "V is a subspace in "H^. 
Suppose that the following conditions Ci — C4 hold: 

Ci . The linear bounded operators R'^ : 'H^ — > 'H'^ exist such that \\R'^f\\^^ — tII/II^o/'"' (^ny 
f e'V. Here y > is a constant. 

Cj- Operators are positive, compact and self-adjoint. The norms ||>^'^||xch'^) are 

bounded uniformly in e. 

C3. For any f e "V: m'^R'^f - R'^Ji^fW^s 0. 

e->0 

C4. For any sequence f^ € 'H^ such that sup < 00 the subsequence e' c e and w e 'V 

£ 

exist such that \\Jl^f^ - R^w\\^s — > 0. 

e=e'^0 

Then for any ^ e N 
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where (ind {/^/t}^i are the eigenvalues of the operators and J^, which are renumbered 

in the increasing order and with account of their multiplicity. 

Let us apply this theorem. We set 'H^ = ^2(6^, -H^ = LjCR) ® LzCBy), = (-Ag, + 

= {£j + l)-\'V = -K". We introduce the operator R"" ■.'H° ^ "K^ by the formula: 



[R'f]{x,z) = 



if\x). 



[/^(x,z-Mcos0p, (x,z)€Bp 



(x,0)€R'^, PR „ 

' / = (A/^)€'K°-L2(R)©L2(By) 



We also denote H^{R) = [u e H\R) : mI^r - o), H^iG^j) - [m € H^Cj) : m|5gj ^ O), -K' = 
//q(R) ® -?^'(Bj) c "K" and introduce the operator Q'^ : H^^iG^j) 'H^ satisfying the properties 
that are similar to those of the operator (see above): 



V£ > 0, Vv € //^(Gp : R'^Q'v = v, \\Q'v\\^^ < C||v||^,(Gp 



(3.30) 



Evidently conditions Ci (with 7=1) and C2 hold. We verify condition C3. Let / € 'K'^. We 

set /'^ = R'^f, = J^l'^/^, = Q^v". One has 

J ((Vv^Vw'=) + M'^w^-/^w'^)^/V = 0, Vw^e//^(Gp (3.31) 



Clearly the norms Hv'^lP , ^ are bounded uniformly in £. Taking into account (3.301) we con- 
elude that the subsequence (still denoted by e) and v = (v^, v^) e TY' exist such that 



£-»0 



V weakly in "K^ and strongly in "K" 



Let w € -K^ = [w - (w^, w^) e "K' : supp/^ c R \ {(0,0)), supp/^ c B^- \ |(0, -h-)^, i.e. 

= in a neighbourhood of {(0, 0)), = in a neighbourhood of {(0, We set w^ = R'^w. 

Then, when e is small enough, = in some neighbourhood of 5B^ and w*^ € //^(Gp. Substi- 
tuting w'^ into (3.31 1 we obtain (e is small enough): 

j ((Vv^^,Vw^) + rV-/V)jx + J ((Vv^^,Vw^) + v^V-/V)jV = (3.32) 



Passing to the limit in (3.32 1 as £ ^ and taking into account that the space "K' is dense in "K^ 
(see e.g. |[25l ). we obtain the equality = v that obviously implies the fulfilment of C3. 

Finally, condition C4 follows from the fact that if sup H/'^H-h^ < 00 then the norms WQ'^'^^f^W'^i 

E 

are bounded uniformly in e and, therefore, the subsequence e' <z s and w € "K' exist such that 
ge^eyfi — ^ ^ weakly in "H^ and strongly in IH^ 

Thus, the eigenvalues /i^ of the operator converge to the eigenvalues p.k of the operator ^ 
as e ^ 0. But /if (Gfp = (pP'^ - 1, /I*: = (jU/t)"^ - 1 that implies (3.29|l. The lemma is proved. 

Lemma 3.4. For j - I,. . . ,m: 



limV . 



CTj-A 



llL2(n) 



(3.33) 
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Proof. For x € G^j we denote F{x) - dhtgc:{x,S'?j), where by distgE(-, •) we denote the distance 
with respect to the metric g^. We introduce the set 

5f.ra = {ye Gf. : fiy) = nx)] 

Obviously S'^j[x] is a (« - l)-dimensional sphere (in particular if ^ e 55?. then F{x) = ks/2 and 
Sf.[x]=5Bf.). 

We define the function m?.(^) by the formula: 

Using the Poincare inequality (for the spheres S ij[x]) we get 

2 



\U — U 



< c2max(diam5f,.ra)2||VMl|^^(C.) < Cs^\\Vu%l^^^.) 



{u^}s'i.- Clearly u^^. e dom(A^s ) and 



(3.34) 



In view of Lemmas 



3.2 



3.3 



and since A ^ Ul'^jl' ^ ^"("^ge) when e is small enough. 

7=1 



Therefore, the foUowmg expansion is valid: 

DO 

^^ = Z4(^)' where /*(£) 



vf(Gf,) 



k=l 



vf(Gf,) 



L2(Gf,) 



(if (Gf .) - i^) 



(3.35) 



Here /^^ = A'^(u'^)s'^., {vf (G^y)}^_^ is a system of the eigenfunctions of -A^^ corresponding to 



We denote A'^ = max max U'^ - /if (G^ ) . Thus, it follows from Lemma 

7=1, m k=2,oo 



Therefore, taking into account (2.7 1 and using Lemma 3.1 we obtain 

2 

DO ^ 



3.3 



that hm A"" ^ 0. 

£^0 



iel-^ k=2 



e" ^ 

£^0 



(3.36) 



As in Lemma 



3.2 



we denote vf. = vf(Gf ). We normalize vf. by the condition {vf^B". = 1- 



Using the estimates! 3.23 1, (3.24i and Lemma 3.2 we obtain that 

2 



llL2(Bf,) 



A'pjM 



|2 



(CT,- - i)2 



(3.37) 



as e — > 0. Thus, it follows from (3.35 1-(3.37 1 that 



lim V uf 



|2 



i2(B^) (cry - i)2 
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(3.38) 



Finally, using (3.35 1, ( 3.36[ ), ( |3.38[ ) and Lemma [3J] we get 

y \k ' - y 



i2(Bf,) 



h 



ix)dV'' + 



(CTj - Ay (T 



lApj 
+ -,+Pj 



(Tj-A 



\B1 



llL2(n) 



£-»() 



(3.39) 



Then (3.33 1 follows from (3.34i and (3.391. The lemma is proved. 
Lemma 3.5. For any w € C^iu) the function e C°°(n) exists such that: 



w + € C"(n) 



max w^'' 







%x)\ + max IVw'^Cjc)! 

Proof. We define the function V e C°°(R") by the following recurrent formulae 

1 (^x\ , . . . , Xji) — A^fii^xi , . . . , Xii—\)Xii + Bf^(^x\ , . . . , Xf^—\\ 

a - 2, . . . ,n : 



(3.40) 
(3.41) 



Ba{Xi,. . . ,Xa-l) - Aq,_i(Xi, . . .,Xa-2)Xa-l + Ba-l{Xi,. . . ,Xa-2), 
Aa{Xi,...,Xa-l) = {0%l9a - \)Ba{Xi, . . . , X^-l), 



Bi = 1, Ai - eyoi - 1. 

It is easy to see that max|l''(;c) - 1| + max|Vl''(x)| and I'' € C~ ,„(□), where := 

xen xen e^O ' 

(^/6li, . . . , Gf,/en). Then we set 

- (1^ - l)w 



Obviously the function satisfies the conditions (3.40 1, (3.41 1. The lemma is proved. 
We continue the proof of Theorem 2.1. For an arbitrary w'^ e dom(77M'!) we have 

j {{Vu^Vw%-A''u^w^)dV' = (3.42) 

where (Vm'^, Vw'^)e is the scalar product of the vectors Vm*^ and Vw'^ with respect to the metric g^. 

We substitute into (3.42) the test function of a special type. Namely, let w be an arbitrary 
function from C~(n), w'^ e C°°{a) be the function satisfying ( |3.40| ), ( |3.41| . Let wj, j ^ \,...,m 
be arbitrary functions from C°°(n). Let 0(r) be a twice continuously differentiable non-negative 
function equal tolasO<r< 1/4 and equal to as r > 1 /2. We set 



\x-x^j\-dy 



0)^. - o 



U-x^-d'^. 

KE 
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Then we set = w"^ + 6"^, where 
w(x), 



^ e \ U UR- 



+ 



(x) + (w(xf.)-w(x))0f.(x) 

+(w/xf^)-w(xpv^(x)cl);^(x), -xeRl 

'y(xf .) + (w(xf .) - W/Xf.)) (l - Vf .(X)) , X € B, 



(3.43) 



(5'=(X) 



w(x), 



X e \ u U ^; 



w 



(X) + (r (xf .) - w\x)) (Df .(X), X € /?f . 

xeS^. 



^^(xf.), 



Herethefunctionvfasdefinedby (3.181, (3.151, (3.16). It follows from (3.40 1 that w'^ € dom(?7M=). 
Substituting this into ( |3.42 i and integrating by parts we obtain 

J {-u^^y|^ -A''u^w'')dV^ + J yK]M^J5^ + J i(Wu\W6% - A''u''5'')dV^ = (3.44) 



where v is the outward normal vector field on ^M'^. 



In view of (2.6 1-(2.7 1 and the Cauchy inequality, the last term in (3.44) is estimated by 



C||m^||^1(Me) ^/max [^(x)!^ + max |Vh'(x)P and tends to zero as £ ^ in view of (3.41 

J 

Jdi 



In view of (|3.8|) the second term tends to j v [w] uds as e ^ 0, where v is the outward 

Ida 

normal vector field on da, ds is the density of the Lebesgue measure on da. 

Now let us investigate the first term. Firstly we study the integrals over Q.^. Integrating by 
parts we get 



m „ 

I -A{(w(xf.)-w(x))Of.(x))M^(x)JV^ 

TP .•_ 1 



re J- j=l 



J (v {(w(xfp - w(x))a)f,.(x)) , V'=n^M'=(x))jx - J AwWu^'dx 



< C(w) ■ WWu' 



in 



(3.45) 



Hereafter by C(w) we denote a constant depending only on w. 
Let us prove that the function e ^2(n)> 



'J 



Z Z {(w^(xf .) - w(xf .))vf .(x)cDf .(X)) , X e /?; 

m 

0, X e □ \ U U ^' 



converges weakly in L2(n) to the function ^ o"7P7(>v - Wy). Indeed using the properties of 
X e R% : Avf,.(x) - 0, |D"vf,.(x)| < Ce"|x - xfJ^^""'"!, a - 0, 1 
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and the enclosure supp(D'^(l)p c |x € : keIA <\x- x^j\ < (a + 0) we obtain 

J |-A {(w/xf.) - w(xf .))vf -WOf .(x))|' dx < C(w)£" (3.46) 

Hence the norms ||^'^||L2(n) are bounded uniformly in e. Taking into account (3.46 1 we obtain for 
an arbitrary / € C^iu) (below is the normal vector field on (9D^^. directed outward 7?^^: 

m p 

T,T, -Aj(w/xfp-w(xpvf.(x)cl)f.(x))/(x)jy- = 

ij 

m „ 

-ZZ-^(4X^(4)-^X4)) v^K]^5-+o(i)- 

'J 

m m „ 

= Z Z /(4)(^(4-) - ^;<4-))c^,-p,e" + o(l) ^^^j 2] (TjPj \ /(x)(w(x) - w/x))^/x (3.47) 

:^ TP. .■_ 1 1 



Here we have used the following computations (below r -\x- x^.|): 



dr 



{dp"-'ajn-i ~ >„-i(n - 2)d'j-^E" = (Tjpjs", E^O (3.48) 



Since C°°(a) = L2(n), then converges weakly in L2(n) to ^ o-jPj(w - wj) as e ^ 0. 
Using this, (3.5 1, (3.6 1 and (3.45 1 we conclude that 

-Aw'^u'^dV'^ = J -Awu + ^ cTjPjiw - Wj)i 



7=1 



dx 



(3.49) 



In the same way (using the estimate ( 3.20 1) one can prove that 

lim I / 



lim^J A^w'^u'^dV'' 
as 



Awudx 



Now, we investigate the behaviour of the integrals in (3.421 over |J B"^.. Using (3.19 



asymptotics), (3.48 1 and the Poincare inequality we get 



'J 



(3.50) 
(the last 



22 r -A[(w(xf.)-w/xf.))(l-vf/x))]«'=(x)JV^ = 

'J 

m i~. 

- Z Z<"'H(^^-^ - -v^K]^5^ + o(i)^ 

'^i'' air: 

IJ 

m p m p 

-'^^o'jPj I [Wy - H'](x)nyM'^(x)(ix + o(l) — Q-jPj I (w/x) - H'(x))My(x)(ix (3.51) 
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where [w /- w] € L2(n) is a step function: [wj-w]{x) = H';(xf.)-w(xf.), x € n^, / € J^; it is clear 

J J J ij ij I 

that [wj - w] converges to wj - w strongly in L2(n) as e ^ 0. 



In a similar manner we obtain 



m „ m „ 

hm ^ ^ I '^''^''"'''^^^ ~ "^^^Pj I Wj^jdx 

^"^iel^j=li ;=1 i 



(3.52) 



Thus, from (3.49l-(3.52) we obtain that the functions u e dom{T]^), uj e L2(n) {j = l,...,m) 
satisfy the equality: 



-AWM + ^ (TjPjUiw - Wj) + ^ (TjPjUj{Wj - W) 
7=1 7=1 



dx + J' V [w] uds- 

dn 



! 



7=1 



dx^O (3.53) 



for arbitrary w e C^(n), wj e C~(n) (7 - 1, ... , m). 

I we c 

-, k = 1, . . . ,m 



Substituting w = 0, = 0, j t k into (3.53 1 we obtain 

(TkU 



Uk 



(Tk-A 



(3.54) 



Then substituting into (3.53 ) wj = (V7), integrating by parts and taking into account (3.54 1, we 



conclude that u e dom(?7Q) satisfies the equality 



J [(Vu,Vw) - Ar{A)uw]dx ^ 0, VweC^in) 



where fiA) is defined by (2.101. Hence u € dom(AQ) and 



-A> - ATiA)u 



In view of Lemma |3.4| m 0. Then A'FiA) e cr(-ARn) and, therefore, due to ( |3.3[ ) A e cr{^) \ 

m 

7=1 



The fulfilment of property (|A|) is proved. 
3.3. Properly of Hausdorjf convergence 



Let A € o-(^) n [0, L], L ^ IJ We have to prove that there exists A"" € o-(-Aa/<:) n [0, L] 

such that A'^ A. 

At first we prove property (|B]) for the case A < L. 

We assume the opposite: the subsequence (still denoted by e) and 6 > exist such that 



dist(/l,cr(-AM'^)) > 6. 



(3.55) 



Since A e <t(JI), then the function F = 



/i 



€ // exists such that 



F ^ im(J?l - AY), where I is the identity operator 
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(3.56) 



Let f^ix) € L2{M'^) be defined by the formula 

[fix), X € Q.\ 



fix) = 



It follows from the Cauchy inequality and (2.7 1 that the norms ||/'^||l2(M'^) are bounded uniformly 
in £. 

Inequahty (|3.55|) implies that /} € R \ o-{-Am'=). Then im(-AM= - /II) = L2{M'^) and thus, the 



unique e dom(AM=) exists satisfying 



(3.57) 



In consequence of (3.55 1 u'^ satisfies the inequality 



\\U%L^(M-) < 6-'\\f\\L.AM^) < C 

Furthermore 

Then there exists a subsequence (still denoted by e) such that 

U^u'' ue H\M.") weakly in H\E.") and strongly in L2(G) for any compact set G c R" 
Wju"" Uj € L2(R") weakly in L2(R") (j - 1, . . . , m) 

where n^, IVj (j = 1, . . . , m) are the extension operators introduced in the previous subsection. 
For an arbitrary function w*^ € C^(M'^) we have 



J' 



(3.58) 



Let w e C^(R"), wj € C[^(R") (J = \,...,m) be arbitrary functions. Using them we construct 
the test-function by formula (3.43 1 (but with R" instead of and with Z" instead of I^) and 



substitute it into (3.58 1. Performing the same calculations as in the previous subsection we obtain 



{Vu, Vw) + ^ 0-jPjU{w - Wj) + ^ (TjPjUj{Wj - w)- 



7=1 



7=1 



-A 



"^ + 2]P7"7^7 - fw + Y^PjfjWj 
7=1 ) \ 7=1 



dx^O (3.59) 



for arbitrary w e C^(R"), wj e C^(R") (7 = 1, . . . , m). It follows from (3.59 1 that 



U = 



yUfi 



€ dom(:?l) and J{U - AU ^ F 



We obtain a contradiction with (3.56 1. Then there is A'^ e cr(-AM''') such that lim/l'' = A. Since 



A < L, then A'^ < L when e is small enough. 



£-»0 



Finally, we verify the fulfilment of property (B 1 for the case A = L. Since L ^ |J ju'; , then 

1 I j^i^ Ji 



(3.1 1 implies that {L - d,L - 6/2) c cr{Jl) whend is small enough. Let As € (L - S,L - 6/2). 
We have just proved that if e < £{6) then A^ e o"(-Am'^) exists such that \A^ - As\ < 6/2. Then 
A'^ € (L- 36/2, L) as e < £{6) that obviously implies the fulfilment of property ([B]). 
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3.4. End of the proof 

In the proof of the HausdorfF convergence we used the fact that M'^ is F-periodic manifold, 
is a period cell. Now let us recall that is also F'^-periodic manifold, is a corresponding 
period cell (/ is arbitrary, so from now on we consider / - 0). Then 



(T{-^M^) = [j[al,bl^ 



where [a^, Z?^] = {/i^(Mp, 9 e T" 
Lemma 3.6. lim b 



ni+1 



Proof. As usual by /I^(Mq) we denote the k-th eigenvalue of the operator -A^^, which is the 
Laplace-Beltrami operator on Mg with Neumann boundary conditions. 

Using the same idea as in the proof of Lemma 3.3 (i.e. e~^-homothetic image of M^), we get 



lime-^ii(M^) = Ak, k = 1,2,3... 



(3.60) 



where {/l/tlieN the eigenvalues of the operator X which acts in the space L2(n) ® ^2(B;) and 

y=l,m 

is defined by the operation 



X = - 



Ab. 



AS 



I 



^ 




Recall that □ is the unit cube in W, By is the n-dimensional sphere of the radius bj (j = 1, . . . , m). 
One has Aj = Ai(Bj) = 0, j = 1, . . . , m, Am+i = A^{n) = 0, and 

A^+2 = min {^^ (□), AjiBj), j = l,...,m]>0 



Thus, in view of (3.60 1 lim A'^^^^^{M^^) = oo. Due to inequality (1.3 1 A'^^^^^iM^^) < a^^^- Thus, 
lim a 



'm+2 



Suppose that there exists a subsequence (still denoted by e) such that the numbers ftj^^^j are 
bounded uniformly in e. Let L > max /u and L > b^^^^^. Let Li > L. Since a^^^^^ — then 



;=l,m 



^m+2 ^ ^1 when e is small enough. Hence o"(-Am=) n [L, LJ = when s is small enough. But 

oo. The lemma is 



this contradicts to property ( B I of the Hausdorff convergence. Hence b 
proved. 



m+l 



It follows from Lemma 3.6 that within an arbitrary finite interval [0, L] the spectrum cri-i^M") 
has at most m gaps when e is small enough, i.e. 



cT{-^M■^) n [o,L] = [o,L] \ \J{(T),^fJ) 



(3.61) 



where (cr^, c [0, L] are some pairwise disjoint intervals, nf < m. Here we renumber the 
intervals in the increasing order. 
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Let L > max ju , be arbitrarily large number. We have just proved that as e ^ the set 

i=\,m 



(tC-A^e) n [0, L] converges to the set cr{Jl) n [0, L] = [0, L] \ ^ |J {crj,^j) j in the Hausdorff sense. 



Then by Proposition JiAnf - m when e is small enough and 

yj=\,...,m: limcr^^cr;, limff-=fij 



Finally, we denote by J''^ the union of the remaining gaps (if any). Since Zj^^^j — > oo and 



b"" , < inf J"", then 

m+l ^ ' 



inf J'^ > L 



when e is small enough. This concludes the proof of Theorem 2.1 



Remark 3.1. Actually, we have proved a slightly strong result: lirna^^^ = ^k, lirn^^ - (^k, 
k = 1, . . . , m, lim , = oo, i.e. the first m gaps of the spectrum ctC-Am*^) (e is small enough) are 
located exactly between the first (m+l) bands. 



4. End of the proof of Theorem 0.1 ; choice of the constants dj, bj and conclusive remarks 



In order to complete the proof of Theorem 0. 1 we have to choose the constants dj, bj in {2.6), 



(2.7 1 such that equalities (0.8 1 hold 



Theorem 4.1. Let {aj,/3j) (j= 1, . . . ,m, m e M) be arbitrary intervals satisfying {0.4). Let M"^ 
(e > 0) be an n-dimensional periodic Riemannian manifolds of the form (|2.i|). 



Then {0.8) holds if we choose 



dj 



ajn-\{n - 2) 

i=l,m|(# j 
i=l,m\i^ j 



LJ \ai-aj 



n>2 



n - 2 



'nOij . J_i \ai- ajj 



i=l,m\ii: j 



i=l,m\i^ J ^ ' J 



(4.1) 



(4.2) 



Remark 4.1. Since the intervals {aj,/ij) satisfy (0.4 1, then 

Vj : aj < Pj, V/ i= j : sign(yS,- - Uj) ^ sign{ai - aj) i= 

Therefore, the expressions (J3j - aj) Yi \ | > j = 1 , . . . , m are positive and thus the 

choice of dj and bj is correct. 



Proof. Substituting dj, bj (4.1 1, (4.2) into (2.8 1 we get 
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i.e. the first equality in (0.8 1 holds. Furthermore substituting bj (4.2 1 into (2.9 1 we obtain 



PJ 



n 

i=l,m|(5t j 



a: - a, 



(4.3) 



It remains to prove that luj = Pj. Recall that the numbers fij (j = I, . . . ,m) are the roots of the 



equation (2. 10 1. Therefore, in order to prove the equality lUj = Pj, we have to show that 



l,...,m: 2 



(4.4) 



Let us consider ( 4.4 1 as the linear algebraic system of m equations with unknowns pj (j 
1, . . . , m). In order to end the proof of theorem we have to prove the following 



Lemma 4.1. The system has the unique solution pi,...,pm which is defined by {4.3). 

Proof. We prove the lemma by induction. For m = \ its validity is obvious. Suppose that we have 
proved it for m = A'^ - 1 and let us prove it for m - N. 

Multiplying the k-th equation in (4.4 1 (^ = I, . . . ,N)hy j3k - and then subtracting the A^-th 
equation from the first N - I equations we obtain 



Mk=\,...,N -\: y "^^^ = 



1 



where the new variables py, 7 = 1 , . . . , - 1 are expressed in terms of pj by the formula 

On - aj 



Pj ■= Pj 



^, j=\,...,N-\ 



(4.5) 



Thus, the numbers pj satisfy the system (4.4 1 with m = N - 1. Therefore, by the induction 

^j-"j n IPi-oij' 



Pj 



a. 



n 



Of - a 



(4.6) 



It follows from (4.5 1, (4.6 1 that pj, j = I, . . . ,N - I, satisfy formula (4.3 1. The vaUdity of this 
formula for p^ follows from the symmetry of the system. 

Lemma 4.1 and Theorem 4.1 are proved. This completes the proof of the main theorem. 

Remark 4.2. We noted above that the metric g'^ of the manifold M'^ is continuous but piecewise- 
smooth (see formulae ( 2.3 1-( 2.4 1). However one can approximate by a smooth metric g^'' which 
difi'ers from only in small p-neighbourhoods of dB'lj and moreover the corresponding Laplace- 
Beltrami operator has the same spectral properties as £ ^ 0. 

Namely, in a small neighbourhood Ufj of dB'^j we introduce the local coordinates (xi, . . . , x„) 
by formulae ( 2.2 1 and define g^'^ by the formula 

Sl^pixi, ...,X„) = gl^^iXl,. . . , Xn)ip{Xnlp) + gt„p{xi, . . . , X„)(l - (fiiXn/p)) 

where ^(r), r € R, is a smooth positive function equal to 1 as r > 1, equal to as r < -1 and 



positive as -1 < r < 1, the coefficients g%^^ are defined by (|2.4|). Outside |J Uf. we set = g^. 
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It is easy to see that A^^g^ < g^f < B'^Pg'^, where A'^p, B'^p are positive constants depending on 
s and p in such a way that for fixed e 



MmA^'P = limB^P = 1 



(4.7) 



Using the min-max principle one can obtain that 



(•gep)l+n/2 



i^(Mf)<4(Mf,/0< 



(B^'P) 



nil 



4(Mf) (4.8) 



Here /l^(Mp g'^P) is the A;-th eigenvalue of the Laplace-Beltrami operator with ^-periodic boundary 
conditions on the manifold equipped with the metric g^P. This inequality is proved in HI 
Chapter A] for manifolds without a boundary, for our case the proof is completely analogous. 
Let 6\ > 0, Li > 0. We have just proved (see Theorems 2.1[ 4.1 1 that there are such e = 



e{6i,L]) and such dj,bj that the manifold M = satisfies (0.2l-(0.3 1 with 6 = 6i, L = Li. 



So let us fix e ^ e(Li, (5i). Then it follows from (4.7 ), (4.8 1 that 



V6» € T", V/t € N : |4(Mf) - 4(Mf,/^)| ^ 

IK! K I I 



(4.9) 



}, where G is any compact subset of N. Then using ( 1.4i, (4.9 1 



uniformly in (9, k) from T" 
and taking into account Remark [XT] we conclude: there is such p - p{e{5\,L\)) that the manifold 
(M^,gep) satisfies (|a2]l-(|03]l with 5 = 25y, L = U - 6^. 



Now, let 5 > 0, L > 0. Setting 5\ = 6/2, Li = L + 6/2 we conclude that the manifold 



M = {M'',g'^P), where s = s{6uU), p = p{s{6uU)), satisfies (0.2i-(0.3l. 



Acknowledgements 

The author is grateful to Prof. E. Khruslov for the helpful discussion. The work is supported by 
the French-Ukrainian grant "PICS 2009-201 1. Mathematical Physics Methods and Applications". 



References 

[1] C. Anne, Exemples de convergence de valeurs propres sur des surfaces ayant una anse tres fine. These, Universite 
S. et M. Grenoble, 1984. 

[2] L. Boutet de Monvel, E.Ya. Khruslov, Averaging of the difi"usion equation on Riemannian manifolds of complex 

microstructure. Trans. Mosc. Mat. Soc. (1997) 137-161. 
[3] L. Boutet de Monvel, I.D. Chueshov, E.Ya. Khruslov, Homogenization of attractors for semilinear parabolic 

equations on manifolds with complicated microstructure, Ann. Mat. Pura Appl. (4) 172 (1997) 297-322. 
[4] Y. Colin de Verdiere, Construction de laplaciens dont une partie finie du spectre est donnee, Ann. Sci. Ec. Norm. 

Super. (4) 20(4) (1987) 599-615. 
[5] G. Dal Maso, R. Gulliver, U. Mosco, Asymptotic spectrum of manifolds of increasing topological type. Preprint 

S.I.S.S.A. 78/200 1/M, Trieste, 2001. 
[6] E. B. Davies, E. M. II Harrell, Conformally flat Riemannian metrics, Schrodinger operators, and semiclassical 

approximation, J. Difi'er. Equ. 66(2) (1987), 165-188. 
[7] P. Exner, O. Post, Convergence of spectra of graph-like thin manifolds, J. Geom. Phys. 54(1) (2005) 77-1 15. 
[8] A. Figotin, P. Kuchment, Band-gap structure of the spectrum of periodic dielectric and acoustic media. I. Scalar 

model, SIAM J. Appl. Math. 56(1) (1996) 68-88. 
[9] A. Figotin, P. Kuchment, Band-gap structure of the spectrum of periodic dielectric and acoustic media. II. Two- 
dimensional photonic crystals, SIAM J. Appl. Math. 56(6) (1996) 1561-1620. 
[10] N. Filonov, Gaps in the spectrum of the Maxwell operator with periodic coefficients, Commun. Math. Phys. 

240(1-2) (2003) 161-170. 

[11] L. Friedlander, On the density of states of periodic media in large coupling limit, Commun. Partial Differ. Equ. 
27(1&2) (2002) 355-380. 



29 



[12] E.L. Green, Spectral theory of Laplace-Beltrami operators with periodic metrics, J. Differ. Equ. 133 (1997) 15-29. 
[13] R. Hempel, I. Herbst, Strong magnetic fields, Dirichlet boundaries, and spectral gaps, Commun. Math. Phys. 
169(2) (1995) 237-259. 

[14] R. Hempel, K. Lienau, Spectral properties of periodic media in the large coupling limit, Commun. Partial Differ. 

Eq. 25(2000) 1445-1470. 

[15] R. Hempel, O. Post, Spectral Gaps for Periodic Elliptic Operators with High Contrast: an Overview, Progress in 

Analysis, Proceedings of the 3rd International ISAAC Congress BerUn 2001, Vol. 1, 577-587, 2003. 
[16] G.A. losifyan, O.A. Olejnik, A.S. Shamaev, On the limiting behaviour of the spectrum of a sequence of operators 

defined on different Hilbert spaces, Russ. Math. Surv. 44(3) (1989) 195-196. 
[17] A. Khrabustovskyi, Asymptotic behaviour of spectrum of Laplace-Beltrami operator on Riemannian manifolds 

with complex microstructure, Appl. Anal. 87(12) (2002) 1357-1372. 
[18] A. Khrabustovskyi, On the spectrum of Riemannian manifolds with attached thin handles, J. Math. Phys. Anal. 

Geom. 5(2) (2009) 145-169. 

[19] A. Khrabustovskyi, Homogcnization of eigenvalue problem for Laplace-Beltrami operator on Riemannian mani- 
fold with comphcated 'bubble-like' microstructure. Math. Methods Appl. Sci. 32(16) (2009) 2123-2137. 

[20] A. Khrabustovskyi, Homogcnization of spectral problem on Riemannian manifold consisting of two domains 
connected by many tubes, submitted for publication; arXiv:1011.3931. 

[21] F. Lledo, O. Post, Existence of spectral gaps, covering manifolds and residually finite groups. Rev. Math. Phys. 
20(2) (2008) 199-231. 

[22] V.A. Marchenko, E.Ya. Khruslov, Homogcnization of Partial Differential Equations, Progress in Mathematical 
Physics 46, Birkhauser, Boston, 2006. 

[23] O. Post, Periodic manifolds, spectral gaps, and eigenvalues in gaps, Ph.D. thesis, Technische Universitat Braun- 
schweig, 2000 (http://www.mathematik.hu-berlin.de/''post/papers/promo-e.pdf). 

[24] O. Post, Periodic manifolds with spectral gaps, J. Differ. Equ. 187(1) (2003) 23-45. 

[25] J. Rauch, M. Taylor, Potential and scattering theory on wildly perturbed domains, J. Funct. Anal. 18 (1975) 27-59. 
[26] M. Reed, B. Simon, Methods of Modem Mathematical Physics IV: Analysis of Operators, Academic Press, 1978. 
[27] L. Tartar, The general theory of homogenization. A personalized introduction. Springer, Berlin, 2009. 
[28] M. Taylor, Partial Differential Equations I. Basic Course, Springer- Verlag, Heidelberg, 1996. 
[29] V.V. Zhikov, S.M. Kozlov, O.A. Oleiruk, Homogenization of Differential Operators and Integral Functionals, 
Springer, New York, 1994. 

[30] V. Zhikov, On spectnmi gaps of some divergent elliptic operators with periodic coefficients, St. Petersb. Math. J., 
16(5) (2005) 773-790. 



30 



